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ABSTRACT 

We investigate the relationship between the dark matter and baryons in the linear regime. 
This relation is quantified by the so-called "filtering scale". We show that a simple gaussian 
ansatz which uses the filtering scale provides a good approximation to the exact solution. 

Subject headings: cosmic microwave background - cosmology: theory - cosmology: large-scale structure 
of universe - galaxies: formation - galaxies: intergalactic medium 

scale becomes essentially irrelevant, because grav- 
itational instability leads only to slow, power-law 
growth of fluctuations. Let us consider the fol- 
lowing simple thought experiment: in a universe 
with linear dark matter fluctuations, the gas is in- 
stantaneously heated to high temperature. The 
Jeans scale also increases by a large factor in- 
stantaneously, whereas it takes about a Hubble 
time for the fluctuations in the gas to respond to 
the changed pressure force. Thus, the instanta- 
neous value of the Jeans scale does not correspond 
to the characteristic scale over which the fluctua- 
tions are suppressed, but instead the suppression 
scale - which we call the "filtering" scale follow- 
ing Gnedin & Hui (1998) - depends on the whole 
previous thermal history. Only in the unphysical 
case of temperature evolving as an exact power- 
law of the scale factor at all times does the Jeans 
scale become proportional to the filtering scale (Bi, 
Borncr, & Chu 1992; Fang et al. 1993). However, 
incorrect expressions for the pressure force filter- 
ing have been used even until quite recently (c.f. 
Choudhury, Padmanabhan, & Srianand 2001) 

In this paper we investigate the role of the fil- 



1. Introduction 

If humans had dark matter vision, cosmology 
would be a solved problem by now. However, this 
is not the case, and in our attempt to understand 
the distribution of matter in the universe we need 
to rely on baryons - stars and gas - that trace the 
underlying distribution of the dark matter. 

The baryons, though, are not the perfect tracer, 
since they are subject to other forces in addition to 
gravity. In this paper we restrict our attention to 
the linear regime, but even in this simplest phys- 
ical situation the relationship between the dark 
matter and baryons (i.e. cosmic gas) is nontriv- 
ial, because gas pressure will erase fluctuations on 
small scales. 

The effect of gas pressure on small fluctuations 
is canonically described by the Jeans scale, which 
is defined as the scale at which the gravity force 
equals the gas pressure force. On large scales 
gravity wins, and small fluctuations grow expo- 
nentially, while on small scales, gas pressure turns 
all fluctuations into sound waves. 

However, in the expanding universe the Jeans 
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tcring scale further. Specifically, it has been sug- 
gested (Gnedin & Hui 1998; Gnedin 1998) that 
the relationship between the gas density fluctua- 
tion 5b and the dark matter fluctuation Sx in the 
Fourier domain can be approximated by the fol- 
lowing simple expression: 

ox 

where k is a wavenumber, and /cf is the wavenum- 
ber corresponding to the filtering scale. Our goal 
will be to investigate the range of validity of this 
approximation. 

2. Linear Dynamics in the Expanding Uni- 
verse 

The linear evolution of fluctuations in the ex- 
panding universe containing dark matter and cos- 
mic gas is described by two second-order differen- 
tial equations (Peebles 1980), 

^ + 2i^ = ^Gp(fxSx + f B 6 B ), 
^ + 2tf^ = 4nGp(f x 5x + f B S B ) 

-&k*6 B , (2) 

where 5x(t,k) and Ssit, k) are Fourier compo- 
nents of density fluctuations in the dark matter 
and cosmic gas, which have respective mass frac- 
tions fx and /s, H(t) is the Hubble constant, a(t) 
is the cosmological scale factor, p(t) is the average 
mass density of the universe, cs(t) is the sound 
speed in the cosmic gas (where the sound speed is 
simply defined by c| = dP/ dp, assuming an equa- 
tion of state that relates the P and p), k is the 
comoving wavenumber and t is the proper time. 

On large scales (k — > 0) the relationship be- 
tween the fluctuations in the gas and in the dark 
matter can be expanded as 

where kp is the filtering scale and in general is a 
function of time. 

The filtering scale is related to the Jeans scale 
kj, 

kj = —^G~p, (4) 



by the following relation: 

1 1 [* , D + (t') + 2H(t')D + (t>) 

k%(t) " D + (t)J kj(t') 




where D + (t) is the linear growing mode in a given 
cosmology (Gnedin & Hui 1998). 

Inspection of equation (5) shows that the filter- 
ing scale as a function of time is related to the 
Jeans scale as a function of time, but at a given 
moment in time those two scales are unrelated and 
can be very different. Thus, given the Jeans scale 
at a specific moment in time, nothing can be said 
about the scale over which the fluctuations in the 
gas are suppressed. It is only when the whole time 
evolution of the Jeans scale up to some moment 
in time is known that the filtering scale at this 
moment can be uniquely defined. 

In order to investigate the properties of equa- 
tions (2) as a function of the thermal history of 
the universe, we parameterize the evolution of the 
sound speed (as a function of cosmological rcdshift 
z) in the following way. Before the beginning of 
reionization at z = zbeg the sound speed is much 
smaller than the sound speed in the photoionized 
intergalactic medium (IGM), and is set to zero. 
Between zbeg and zrei the temperature in the 
IGM rises linearly with redshift. The sound speed 
reaches a maximum at the moment of reionization 
-Zrei with the value of Cs ; rei, after which moment 
it falls off as (1 + z) fi l 2 . This parametrization is 
schematically illustrated in Figure 1. Thus, our 
thermal history is parametrized with 4 parame- 
ters: zbeg, zrei, CS3EI, which we replace with 
the mean temperature at reionization in units of 
10 4 K, T 4j rei, and (3. We adopt a flat cosmolog- 
ical model as our background cosmology, which 
introduces two more parameters: f2o and £l B - We 
adopt the following values as our fiducial model: 



n 


= 0.3 


n B 


= 


Tk.REI 


= 2.5 


ZREI 


= 7 


zbeg 


= 10 


(3 


= 1. 



(6) 



Solutions to equations (2) at three different red- 
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Fig. 1. — Schematic representation of our 
parametrization of the thermal history. 



shifts are shown in Figure 2 with thin lines. As 
can be expected, on large scales (k — > 0) the gas 
follows the dark matter, whereas on small scales 
oscillations in the gas turn into slowly decaying 
sound waves. We also show with thick lines the 
gaussian approximation (1) to illustrate its level of 
accuracy for our fiducial model. Clearly, the gaus- 
sian approximation does not reproduce the small 
scale oscillations, but it appears to do a decent 
job in reproducing the transitional region between 
large and small scales. 

3. Results 

In this section we present our investigation of 
the accuracy of the gaussian approximation (1). 
Our approach to testing the approximation is mo- 
tivated by its applicability. We can envision using 
the approximation as a simpler way of computing 
the rms density fluctuations in the gas: 




10 10 2 

k (h Mpc" 1 ) 

Fig. 2. — Solutions to equations (2) for our fiducial 
model at z = 4 (light grey), z = 1.5 (medium 
grey), and z — (dark grey): thin lines show the 
exact solutions, thick lines give the approximation 
(!)• 
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S B (t,k) 
S x (t,k) 



k 2 dk 



* i r p x^ e ~ 2k2/k2At)k2dk > ( ? ) 

where Px{k) is the dark matter power spectrum. 
Another use of the approximation could be in giv- 
ing the shape of the gas power spectrum on small 
scales: 



P B (k,t) = P x (k) 



6 B (t,k) 



Px(k)( 



5x{t,k) 
2k 2 /k 2 F (k) 



(8) 



Thus, our goal is to estimate the accuracy of the 
gaussian approximation both in equation (7) and 
in equation (8) in our six-dimensional parameter 
space. 

A full sampling of this parameter space would 
be unrealistic, so we focus on varying two of the 
parameters at a time, keeping the rest of the pa- 
rameters fixed to our fiducial values. For the 
dark matter power spectrum Px{k), we adopt a 
canonical Cold Dark Matter + Cosmological Con- 
stant power spectrum with Q^li — 0.25, Harrison- 
Zcl'dovich scalc-frcc primordial spectrum, and the 
BBKS transfer function (Bardeen et al. 1986). 

The choice of the power spectrum does mat- 
ter. For example, if the power spectrum contains 
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Fig. 3. — The accuracy of the gaussian approx- 
imation in computing the rms gas fluctuation as 
a function of zrei and zbeg at z = 4 (grey) and 
z = (black). 
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Fig. 5. — The same as Fig. 3, except as a function 
of fix and zrei (with zbeg begin kept fixed at 
10). 




Fig. 4. — The same as Fig. 3, except as a function 
of T 4>R ei and (3. 



Fig. 6. — The same as Fig. 3, except as a function 
of Qx and 
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Fig. 7. — The same as Fig. 3, except as a function 
of zrei and f3. 

a lot of power on very small scales (below the fil- 
tering scale) , then oscillations shown in Fig. 2 will 
be amplified, which, in turn, will compromise the 
approximation (1). However, we believe that the 
dark matter power spectrum is known sufficiently 
well up to scales of interest, so our choice is justi- 
fied. 

Figures 3-6 show the accuracy of the gaussian 
approximation as a function of our parameters. 
One can see that in general the approximation 
holds very well, with the approximate <jb being 
within 10% of the exact value. The approximation 
breaks down at z > 4 in the limit when reioniza- 
tion is quite fast (zbeg ^ 1-2zrei) and the red- 
shift of rcionization is low (zrei < 6). It also gets 
somewhat worse for very low values of f^o- This 
is not surprising, since in those two regimes the 
oscillations at large k are substantial and make a 
non-negligible contribution to the integral (7). 

We have also found that in essentially all cases 
the gaussian approximation gives a better than 
20% (10% in amplitude) fit to the gas power spec- 
trum (8) for k < 0M F - 



k < 0.9kp. It is important to underscore that 
there is no known physical reason why this approx- 
imation works so well, so it should be considered 
as a mathematical coincidence. Notwithstanding, 
the gaussian approximation can be used in semi- 
analytical models of the Lyman-alpha forest and 
early universe, when a full solution to equations 
(2) is not practical. 
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4. Conclusions 



We have shown that the gaussian approxima- 
tion (1) provides a reasonably good fit (better than 
10%) to the rms density fluctuations in the gas, 
and for the shape of the gas power spectrum for 
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